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1. Introduction 

Top-quark production in hadronic collisions has been one of the most studied signal in 
the last twenty years. Up to recent times, ti pair production has been the only observed 
top-quark source at the Tevatron collider, due to its large, QCD-dominated, cross section. 
Processes where only one top quark appears in the final state are known in literature as 
single-top processes. Their cross sections are smaller than the ti pair one, due to their 
weak nature. This fact, together with the presence of large W + jet and ti backgrounds, 
makes the single-top observation very challenging, so that this signal has been observed 
only recently by the CDF [1] and DO [2] collaborations. 
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In spite of its relative small cross section, single-top production is an important sig- 
nal for several reasons (see also refs. [3, 4] and references therein). Within the Stan- 
dard Model, the single-top signal allows a direct measurement of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element V t j, [5] and of the b parton density. Furthermore, the 
V-A structure of weak interactions can be directly probed in these processes, since the top 
quark decays before hadronizing, and its polarization can be directly observed in the an- 
gular correlations of its decay products [6, 7]. Finally, single-top processes are expected to 
be sensitive to several kinds of new physics effects and, in some cases, are the best channels 
to observe them [8, 9, 10]. For all the above reasons, single-top is an important Standard 
Model processes to be studied at the LHC, where the statistics limitation due to the small 
cross section is less severe and differential distributions can also be studied. 

In order to include a reliable description of both short- and long-distance effects into 
the simulation of hadronic processes, it is important to consistently match fixed order 
results with parton showers. Radiative corrections for single-top production have been 
known for years [4, 11, 12, 13, 14, 15, 16, 17, 18], while the implementation of these results 
into a next-to-leading-order Shower Monte Carlo (SMC), namely MC@NL0 [19, 20], is more 
recent [21, 22]. 

In this work we present a next-to-leading order (NLO) calculation of s- and i-channel 
single-top production, interfaced to Shower Monte Carlo programs, according to the POWHEG 
method. This method was first suggested in ref. [23], and was described in great detail in 
ref. [24]. Until now, the POWHEG method has been applied to ZZ pair hadroproduction [25], 
heavy-flavour production [26], e + e~ annihilation into hadrons [27] and into top pairs [28], 
Drell-Yan vector boson production [29, 30], W' production [31], Higgs boson production 
via gluon fusion [32, 33] and Higgs boson production associated with a vector boson (Higgs- 
strahlung) [33]. Unlike the MC0NL0 implementation, POWHEG produces events with positive 
(constant) weight, and, furthermore, does not depend on the subsequent Shower Monte 
Carlo program. It can be easily interfaced to any modern shower generator and, in fact, it 
has been interfaced to HERWIG [34, 35] and PYTHIA [36] in refs. [25, 26, 29, 32]. 

Single top production is the first POWHEG implementation of a process that has both 
initial- and final-state singularities, and so the present work can serve as an example of 
how to deal with this problem in POWHEG. 

This paper is organized as follows. In sec. 2 we collect the next-to- leading order cross 
section formulae and describe the kinematics and the structure of the singularities. In sec. 3 
we discuss the POWHEG implementation and how we have included the generation of top- 
decay products. In sec. 4 we show our results for several kinematic variables. Most of this 
phenomenological section is devoted to study the comparison of our results with those of 
MC0NL0. We find fair agreement for almost all the distributions and give some explanations 
about the differences we found. Some comparisons are carried out also with respect to 
PYTHIA 6.4, showing that some distributions are strongly affected by the inclusion of NLO 
effects. Top-decay effects are also discussed. Finally, in sec. 5, we give our conclusions. 
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2. Description of the calculation 



In this section we present some technical details of the calculation, including the kinematic 
notation we are going to use throughout the paper, the relevant differential cross sections 
up to next-to-leading order in the strong coupling a s and the subtraction formalism we 
have used to regularize initial- and final-state singularities. In this paper, we always refer to 
top-quark production, since anti-top production is obtained simply by charge conjugation. 

Single-top production processes are usually divided into three classes, depending on 
the virtuality of the W boson involved at the leading order: 

1. Quark- antiquark annihilation processes, such as 

u + d^t + b, (2.1) 
are called s-channel processes since the VF-boson virtuality is timelike. 

2. Processes where the top quark is produced with an exchange of a spacelike W boson, 
such as 

b + u^t + d, (2.2) 

are called t-channel processes. 

3. Processes in which the top quark is produced in association with a real W boson, 
such as 

b + g^t + W. (2.3) 

These Wt processes have a negligible cross section at the Tevatron, while at the LHC 
their impact is phenomenologically relevant. The calculation of NLO corrections to 
Wt processes is also interesting from the theoretical point of view, since the definition 
of real corrections is not unambiguous [22] . 

In this paper we consider only s- and ^-channel processes. In these cases, the POWHEG 
implementation needs to deal with both initial- and final-state singularities, and is thus 
more involved than in processes previously considered. The associated Wt production has 
only initial-state singularities and is thus analogous to previous POWHEG implementations. 
We leave it to a future work. 

In the calculation, all quark masses have been set to zero (except, of course, the top- 
quark mass) and the full Cabibbo-Kobayashi-Maskawa (CKM) matrix has been taken into 
account. However, for sake of illustration, we set the CKM matrix equal to the identity in 
this section. 

We refer to ref. [24] for the notation and for a deeper description of the POWHEG method. 
Here we just recall that with B, V, 1Z and Q we indicate the Born, virtual, real and 
collinear contributions respectively, divided by the corresponding flux factor. The same 
letters, capitalized, are used for quantities multiplied by the luminosity factor. The explicit 
formulae for these quantities are collected in sec. 2.3. 
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2.1 Contributing subprocesses 



In the following, we organize and label all the Born and real subprocesses, keeping the 
distinction between the s- and t-channel contributions. This distinction holds also when 
real corrections are considered, since, due to color flow, interferences do not arise between 
real corrections to s- and t-channel Born processes. 

1) In the s-channel case, there are only Born processes of the type qq' — > tb , where q and 
q' run over all possible different quark and antiquark flavours compatible with the final 
state. We denote with B qq i the (summed and averaged) squared amplitude, divided 
by the flux factor. The corresponding real correction contributions include processes 
with an outgoing or an incoming gluon, i.e. processes of type qq' — > tbg, gq — > tbq' and 
qg — > tbq'. We denote these contributions with lZ qq i, 7Z gq ,( s ) an d 7Z qg ,( s )- Summarizing, 
we have 



process 


notation 


contributing subprocesses 


qq' tb 


Bqq> 


ud — > tb, du — > tb 


qq' tbg 


fcqq' 


ud —7- tbg, du — > tbg 


gq tbq' 




gu — > tbd, gd — > tbu 


qg tbq' 




ug — > tbd, dg — > tbu 



where u and d stand for a generic up- or down-type light quark. 

2) In the ^-channel case, there are only Born processes of the type qb —> tq' (and bq — > tq'), 
where q and q' run over all possible flavours and anti-flavours. Their contributions are 
denoted B q b (Bbq)- We use this notation since we want to keep track of the down-type 
quark connected to the top quark. The structure of real corrections is more complex 
in this case. Contributions obtained from the previous processes by simply adding an 
outgoing gluon, qb — > tq'g, will be denoted as TZ q b- The subprocesses generated by an 
initial-state gluon splitting into a quark-antiquark pair are designated by Tt qg ,(t) for 
qg — > tq'b (lZ gq ^t) f° r 9Q ~ * tq'l>) and lZ g b for gb — > tqq' (lZb g for bg — > tqq'). In the 
former case q and q' are connected via a Wqq' vertex, while the gluon splits into a bb 
pair, so the top quark is color connected with the incoming gluon. In the latter case 
the situation is opposite, since the gluon splits into a qq pair, while the incoming b is 
directly CKM-connected to the top quark. This gives rise to a different singularities 
structure, which we take into account in dealing with the qg —> tq'b (gq — > tq'b) and 
gb — > tqq' (bg — > tqq') processes separately. Summarizing 
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process 


notation 


contributing subprocesses 


bq — > tq' 


B hq 


bu —7- td, bd —> tu 


qb — > tq' 


B qb 


ub —7- td, db —7- tu 


bq -» tq'g 




bu —7- tdg, bd — > tug 


qb tq'g 




ub —> tdg, db — > tug 


gq tq'b 


n gq,(t) 


gu —> tdb, gd —> tub 


qg tq'b 


n qg,(t) 


ug —7- tdb, dg — > tub 


gb tqq' 




gb — > tud 


bg -)• tqq' 


Tlbg 


bg — > tud 



where u and d stand for a generic up- or down-type light quark. 

In order to distinguish s- and t-channel real processes with the same flavour structure, 
we have used the subscript (s) and (t) on the lZ gq and TZ qg contributions. As already 
stated, these contributions do not interfere owing to different color structures, so we can 
keep them distinct. We have drawn a sample of Feynman diagrams for s- and t-channel 
gu — > tdb scattering in fig. 1. 




(a) s-channel (b) t-channel 

Figure 1: Graphs corresponding to s- and t-channel contributions to the real scattering gu — > tdb. 



2.2 Kinematics and singularities structure 
2.2.1 Born kinematics 

Following the notation of ref . [24] , we denote with k Q and k e the incoming quark momenta, 
aligned along the plus and minus direction of the z axis, by k\ the outgoing top-quark 
momentum and by k,2 the other outgoing light-parton momentum. The final-state top- 
quark virtuality will be denoted by M 2 , so that kf = M 2 . The top quark on-shell condition 
is M 2 = m 2 , where nit is the top-quark mass. If K B and K e are the momenta of the 
incoming hadrons, then we have 

fc© = %K @ , (2.4) 
where the momentum fractions, and momentum conservation reads 

k & + k e = h + k 2 . (2.5) 
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We introduce the variables 



and 0±, the angle between the outgoing top quark and the k B momentum, in the partonic 
center-of-mass (CM) frame. We denote with (pi the azimuthal angle of the outgoing top 
quark in the same reference frame. Since the differential cross sections do not depend on 
the overall azimuthal orientation of the outgoing partons, we set this angle to zero. At the 
end of the generation of an event, we perform a uniform, random azimuthal rotation of 
the whole event, in order to cover the whole final-state phase space. The set of variables 
$2 = {s, Y, Q~i, 0i} fully parametrizes the Born kinematics. From them, we can reconstruct 
the momentum fractions 

^ Y - I & -Y 



Xe= V5 e ' Xq= \S 6 ' ^'^ 

where S = {K m + K e ) 2 is the squared CM energy of the hadronic collider. The outgoing 
momenta are first reconstructed in their longitudinal rest frame, where Y = 0. In this 
frame, their energies are 



The two spatial momenta are obviously opposite and both have modulus equal to fc^l^o- 
We fix the top-quark momentum to form an angle 9\ with the e direction and to have zero 
azimuth (i.e. it lies in the xz plane and has positive x component). Both k\ and k 2 are 
then boosted back in the laboratory frame, with boost rapidity Y. The Born phase space, 
in terms of these variables, can be written as 

^4 J c4/^. , u T, I, A <* 3 fc 2 



d<& 2 = dx & dx e (2Tr) 5 (fc e + k e - k\ - k 2 ) 



(2vr) 3 2fc° (2tt) 3 2^ 



where 



= ±JL d sdYdcose 1 ^, (2.9) 

M 2 

(3 = 1 . (2.10) 

s 

We generate the top quark with virtuality M 2 and decay it with a method analogous to 
the one adopted in ref. [37], that will be described in sec. 3.3. We take into account the 
top finite width by first introducing a trivial integration J dM 2 5(M 2 — mf) in eq. (2.9) 
and then by performing the replacement 

5 ( M 2 - mf) ^ . (2.11) 

V ' ^ (M 2 - m 2 f + {m 2 T 2 ) 

With this substitution, the final expression for the Born phase space reads 

d® 2 = l-JL- dM 2 dsdY dcosh^P- ■ (2.12) 
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2.2.2 Real-emission kinematics 

Real-emission processes have an additional final-state parton that can be emitted from 
an incoming leg only {Tl gqM , K gg ,( s ), ^ sg ,(t), fcqg,(t), K gb , lZ bg ) or from both an initial- 
and final-state leg (lZ qq >, TZ b q, Ttqb)- We need then to use two different parametrizations 
of the real phase space, one to deal with initial-state singularities and one for final-state 
ones. We treat the radiation kinematics according to the variant of the Frixione, Kunszt 
and Signer (FKS) subtraction scheme [38, 39] illustrated in ref. [24]. Before giving all the 
technical details, we summarize briefly how the procedure works: 

- We split each real squared amplitude into contributions that have at most one 
collinear (and/or one soft) singularity. 

- We build the collinear (and soft) subtraction terms needed to deal with that singu- 
larity. 

- We choose the emission phase-space parametrization suited for the singularity we 
integrate on. 

In the FKS method, the singular regions associated with the © and e legs are treated with 
the same kinematics. Nevertheless, we have decided to split these two different contribu- 
tions in order to gain a clear subtraction structure. 

We now describe the procedure used to split real squared amplitudes and the corre- 
sponding phase-space parametrizations. Subtraction terms will then be given in sec. 2.3.3. 
We proceed as follows: 

1. We start by considering real processes that have both initial- and final-state emissions, 
namely the TZ qq i, lZ bq and lZ qb contributions. In this case, the FKS parton is the 
outgoing gluon and we choose it to be the last particle. We denote its momentum by 
&3, so that momentum conservation reads 

k S) + k e = k 1 + k 2 + k 3 , (2.13) 

where fe e , k e , k\ and ki label the same particles of the underlying Born process. The 
FKS parton can become collinear to one of the incoming legs or to the other massless 
final-state leg, so we need to introduce a set of functions to project out these different 
singular regions. The general properties these functions have to satisfy were given in 
sec. 2.4 of ref. [24]. In this paper we use 

5 3,e = v -i J_ , s 3 ' e =V- 1 —, S 3 ' 2 =V- 1 —, (2.14) 

«3,e "3,e "3,2 

where ^ 

V = h V - — and dij = h ■ kj . (2.15) 

"3,e «3,e "3,2 

For any given kinematic configuration, these functions satisfy 

S 3,s> + S 3,e + S 3 ' 2 = 1. (2.16) 
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The separation among different singular regions is performed multiplying each real 
contribution with the corresponding S function. For example, for the s-channel lZ qq ' 
case, we have 



-r>3,e 
'\q' 


= Tlqq' 


S 3,e 


<r>3,e 


= 7t qq i 


5 3,e 


7? 3 ' 2 


= Ttqqi 


s 3 > 2 . 



(2.17) 

These contributions are now singular only when the FKS parton becomes collinear 
to k e , k e and k 2 respectively, or soft. Analogous relations hold for TZb q and TZ q b- 

2. Next we consider the real process gb — > tqq' . It is singular when q or q' become 
collinear to the incoming gluon, so that the FKS parton can be respectively q or q' 
and we need again a set of functions to project out the different singular regions. 
Recalling the labeling of the momenta 

g (fe ) b (fc e ) -> t (h) q (k 2 ) q' (k 3 ) , 

we introduce the projecting functions 

S 2 ' m = 



1 


1 ~ 

rf 3,e, 


)" 1 






1 


1 ~ 

^3,e, 


Y 1 1 


d2,s 


/ d 3 , 



5 3 -® = + (2.18) 



to isolate the region where k 2 • k @ — > or ^3 • — > 0. We have then the two 
contributions 

n 3 g f = n gb 5 3 -®, 

K 2 £ = Kgb S 2 >° , (2.19) 

coming from 1Z g b- For bg — > tqq', analogous contributions can be obtained from 
eqs. (2.18) and (2.19) with the substitutions lZ g b — > lZb g and © — > e. 

3. To deal with the remaining real contributions we do not need to introduce any other 
S function, since each of them is singular in one region only (the © one for TZ gq ^ 
and K gq ,(t), the © one for U qg ^ s ) and U qg ^). 

Having split all real contributions in such a way that each term has at most one singularity, 
we can associate with each of them a particular phase-space parametrization, suitable 
to handle that singularity structure. In the following we summarize the reconstruction 
procedure needed to build the real-emission kinematics, given the underlying Born one, 
and a set of three radiation variables. For all the details, we refer to sec. 5 of ref. [24]. 
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Parametrization of the initial-state radiation (ISR) phase space 

The FKS method uses the same phase-space parametrization for describing both the © 
and © singular regions. The set of radiation variables 

•53 = tt. V. # , (2.20) 

together with the Born ones, completely reconstruct the real-event kinematics: $3 = 
(s, Y, 9i, £, y, (/)}. Using eq. (2.7), we can compute the underlying Born momentum frac- 
tions x @ and, from them, we obtain 



9 ~ vr^ V 2 - ^ + y) ' V 2 - m - y) ' [ } 

with the kinematics constraints 

0<£<£ M (y), (2.22) 

where 



£ M (y) = 1 - max 



2(1 + y) 4 



V(l + 5* )"(1 - yf + 16 y a* + (1 - y)(l - *|) ' 
2(1 -y) 4 



(2.23) 



V(l + x|) 2 (! + V) 2 - 16y x| + (1 + y)(l - s|) 
In the laboratory frame, the incoming momenta are given by 

= x & K @ . (2.24) 

In the partonic center-of-mass frame, we define the FKS parton to have momentum 

k' 3 = k' 3 ° (1, sin#sin<^>, sin cos 0, cos 9), (2.25) 

where 

E° = - 

2 

and 



*3° = ^, cos9 = y, (2.26) 



s = (fc e + fe e ) 2 = — . (2.27) 

From eqs. (2.25) and (2.26), we see that the soft limit is approached when £ — > 0, while 
the collinear limits are characterized by y — > 1 (&3 parallel to the © direction) or y — > — 1 
(&3 parallel to the © direction). 

Boosting k' 3 back in the laboratory frame with longitudinal velocity (x® — x e )/ (x e +x e ) 
we obtain k% . Having computed k% and fe® , we can construct k tot = k m + k e — A; 3 , while from 
the underlying Born momenta we have fctot = &i + &2 • We construct then the longitudinal 
boost Ml, with boost velocity /3l = (0,0, /3l), where 

0L = -^i^-, (2-28) 

Xcr, + Xp, 
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so that the boosted momentum k' t ' ot = B^fctot has zero longitudinal component. In addition 
we define 

Pt = -§± (2.29) 

and the corresponding (transverse) boost By, so that ~M>Tk" ot has zero transverse momen- 
tum. The final-state momenta k\ and hi in the laboratory frame are obtained with the 
following boost sequence 

ki = Ml 1 M T 1 M L k i , i = l,2. (2.30) 

Finally, the three-body phase space can be written, in a factorized form, in terms of the 
Born and radiation phase space 

d* 3 = d X@ d Xe {2,)H\k 9 + k e - kl -k 2 - ^) (2 ^ fc (2 ^0 (2 ^ fc = <**£5 , 

(2-31) 

d^ = j^^dtdydcf>=J^{& 2 ,$ l Z) dtdyd<P, (2.32) 



where 



that defines the Jacobian of the change of variables. 

Parametrization of the final-state radiation (FSR) phase space 

For the FSR phase-space parametrization ^ffj 1 , we use the same notation as for the initial- 
state case <J>^|[ (see eq. (2.20)). We define, in the partonic center-of-mass frame, 

7^' y = t^f' <t> = <t>{Vxk, hxkj, (2.33) 

where 

q = k e + k e , k = k 2 + k 3 , (2.34) 

and the notation p stands for \p\ . We denote with ff an arbitrary direction that serves as 
origin for the azimuthal angle of k% around k, while "x" is the cross vector product. The 
notation 4>(vi,v 2 ) indicates the angle between v\ and v 2 , so that <p is the azimuth of the 
vector &3 around the direction k* 

From eq. (2.33) we see that the soft limit is approached when £ — > 0, while the collinear 
limit is characterized by y — > 1 (k^ parallel to k 2 ). 

Given the set of variables <&3 = [s, Y, #i, £, y, <f\ we can reconstruct the full real-event 
kinematics. The momentum fractions x@ are the same as the underlying Born ones, since 
the emission from a final-state leg does not affect them, so that 

= , Xq = Xq and s = s . (2.35) 



*The FKS variant that we use (see ref. [24]) has a slightly different definition of cj> than the one introduced 
in the original FKS papers. 
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Inverting the first identity in eq. (2.33), we immediately have 



where £ is limited by 



with 



k°3=k 3 = ZY> (2 - 36) 



q 2 — M 2 

< C < £m = 5-^ , (2-37) 



M r 2 ec = (g-fc 2 ) 2 = fe2. (2.38) 



The energy (and the modulus) of the other light outgoing parton, always in the partonic 
center-of-mass frame, is given by 

_ q 2 - M 2 CC - 2q% 

fc2 "- 2 " 2 [q o-k 3 (i-y)] ■ (2 ' 39) 

Given k 2 and k% we construct the corresponding vectors fc 2 and ks such that their vector 
sum k is parallel to /c 2 and the azimuth of £3 relative to k (the given reference direction) 
is (p. Having fully defined fc 2 and fc 3 , we can reconstruct the vector k of eq. (2.34) and find 
kmc = q — k. Finally, k\ can be obtained boosting k\ along the k TCC direction with boost 
velocity 

jj — — ( _ — itis£ — rec ^ ^ ^ rec ^2 40) 

\q 2 + (k® cc + k_ Tec ) 2 J k rec 

or, alternatively, exploiting momentum conservation of eq. (2.13). To obtain the momenta 
in the laboratory frame we need to boost back all the outgoing momenta computed in the 
center-of-mass frame. 

In this case too, the three-body phase space can be written in a factorized form in 
terms of the Born and radiation phase space 

= dx B efee^mjfce + k e - kl - k 2 - = d^J , 



(2.41) 



where 



< ;: _ _ 7 s L.-2 f , k 



2 \ -1 



^ d = (4^g^-2?J W ^ 
s 4£ / s£ 



(4vr)3 [2 - £ (1 - y)] 2 V s - M 2 



di dy d<f> = J r F a s d R (* 2) $ r F f d R ) d£ dy d4> . 

'rec / 

(2.42) 

2.3 Squared amplitudes 

In order to apply the POWHEG method, we need the Born, real and soft-virtual contributions 
to the differential cross section, i.e. the squared amplitudes, summed (averaged) over colors 
and helicities of the outgoing (incoming) partons, and multiplied by the appropriate flux 
factor. We have taken the Born, real and soft-virtual contributions from the MC@NL0 code, 



- 11 - 



testing, where possible, our implementation against MadGraph subroutines [40]. All the 
matrix elements have been evaluated in the zero-width approximation, i.e. Ft and Tyy are 
set equal to zero in all the propagators. As already mentioned, to recover finite-width 
effects in top-decay, the top mass M is generated according to a Breit-Wigner distribution, 
centered in m t and with width T t (see eq. (2.11)). 

In the following, we give explicit expressions for the Born and collinear remnant contri- 
butions. Real and soft-virtual matrix elements are more complicated, and we do not report 
them explicitly. Nevertheless, we give the soft and collinear limits of the real amplitude, 
since these expressions are needed in the FKS subtraction formalism. 

2.3.1 Born and virtual contributions 

We denote the s-channel squared matrix element for the lowest-order contribution, averaged 
over color and helicities of the incoming particles, and multiplied by the flux factor l/(2s), 
as B qq i. For example, for the ud — > tb subprocess, we have 



B 



ud 



2s 4 



u(u - M 1 



1 



m 



w 



\V ud \ 2 \V tb \ 2 , 



(2.43) 



where u = (k^ — fo) 2 is the usual Mandelstam variable, g is the weak coupling (e = g sin 9^) 
and Vij's are the CKM matrix elements. Crossing eq. (2.43) we have, for the du initiated 
process, 

1 



M 



m 



W 



\V u d\ 2 \Vtb\ 2 ■, 



and for the t-channel contributions (Bb q and B q b) of the bu — > td and ub 



(2.44) 
td subprocesses 
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\V u d\ 2 \Vtb\ 2 , 



\Vud\ 2 \Vtb\ 2 , 



(2.45) 



where t = (fc e — k\) 2 . The corresponding expressions for bd 
obtained from the latter again by crossing. They are given by 



tu and db — > tu can be 



B h 



1 ° 4 -(- 

2iT n(u " 


M 2 ) 
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t-m 2 w 


i <?v 

— —tit - 
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M 2 ) 


1 





\Vud\ 2 \Vtb\ 2 , 



\V u d\ 2 \Vtb\ 2 . 



(2.46) 



The finite soft-virtual contributions, obtained according to the FKS method, have been 
taken from the MC@NL0 code. We included them in our NLO calculation and tested the 
correct behaviour of our program by comparing our NLO results with the MCFM code [41], 
both for the full NLO cross section and for typical differential distributions. Some com- 
parisons have also been carried out with the program ZTOP [42]. 
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2.3.2 Collinear remnants 

The collinear remnants are given in eq. (2.102) of ref. [24]. Here we limit ourselves to list 
all the contributions, giving only a couple of explicit examples to clarify the notation. 
For the s-channel processes, the collinear remnants are 

ef'(*2, @ ), 0r(*2,e) and 0f(* 2> e), (2-47) 
where the $2,e notation, according to ref. [24], represents the set of variables 

^2,(b = {X( B iX e ,z,k 1 ,k 2 }, with zx^K® + x e K e = ki + k 2 . (2.48) 
The underlying Born configuration $2, associated with the $2,9 kinematics, is defined by 

k 9 = zx m K 9 , k e = x e K e , h = h, k 2 = k 2 . (2.49) 

Similar formulae hold for $2,9- Among the contributions listed in (2.47), only the real 
process qq' — > tbg is singular in both the © and the © region. It thus needs the two collinear 
remnants 



Sf (*2 )@ ) = ^ F (1 + , 2 ) 



1 ^ log^V 10g(1 - z) 



1 — Z J + Z[lj \ 1 — z 
+ (l-z)j B^(l,y,ei). (2.50) 

For the t-channel processes, the collinear remnants are 

0S(*2, e ), 0j(*2, e ), ar(*2,e), SS fl (*2,e), ^(*2,e) and S*(# 2>e ) . 

(2.51) 

In this case, £? ffi fe (3?2,e) contains two terms, since in the scattering gb — > tqq' both the two 
outgoing massless partons q and q' can become collinear to the incoming gluon. We have 

#<•*.) = d - .> d - * + *>) [(^) + . og £ + 2 (!^) J 

+ 2z (1 - z) j [%,(s,YA) +B, 6 (s,y,ei)] , (2.52) 

where Bqn, and £>gft are the corresponding underlying Born processes. All the other contri- 
butions can be obtained in a similar way. 

2.3.3 Soft and collinear limits of the real contributions 

In the FKS formalism, phase-space singular regions are approached when the radiation 
variables £ — > and/or y — >■ ±1. The corresponding singularities are subtracted from 
the real cross section using the plus distributions. One needs to express the singular 
limits in terms of suitable radiation variables and of the corresponding underlying Born 
contributions. In this section we compute these limits and give explicitly their expressions. 
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We start by considering the singular limits of the processes that have both ISR and 
FSR singularities, namely TZ qq >, 1Zb q and lZ q b- These processes are the most subtle, being 
both soft and collinear divergent for initial- and final-state radiation. As an example, we 
study the limits for the s-channel scattering qq' — > tbg. We can deal with ISR and FSR 
separately, having defined the contributions 7t qq ®, T^ qq i and T^ qq i- 

For ISR singularities, we use the set $£ad ^° parametrize the kinematics. When y — > ±1, 
the momentum k% is aligned along the © direction and k% = in the CM frame. The 
real squared amplitude factorizes and we have 



n 



3,® 
-qq' 



47TQq 



P qq (z)B qql = C F 



167ra.q 



(1 + z 2 ) B qq , , 



(2.53) 



where z = (1 — £), P qq (z) is the usual Altarelli-Parisi (AP) splitting kernel and we have 
included the real flux factor l/(2s) and a 1/z factor into the B term, as its definition 
requires. In the FKS approach, one needs the finite quantity £ 2 (1 =F y)T^- q ' q f to perform the 
subtraction of the singularities. In the collinear limit, we have 



3,<i 
qq 1 



y=±l 



C F ^{l + z 2 )B qq , 



(2.54) 



In the same limit, we also note that the contributions [£ 2 (1 =F y)K qq f] and [£ 2 (1 T y)fc qq >] 



go to zero, because the factor ^ 2 (l=Fy) makes them finite and the corresponding S functions 
were chosen to vanish in this limit. 

In the FSR case, the collinear limit is reached when y — > 1. The outgoing momenta kz 
and k 2 become parallel and aligned along their sum, denoted by k. Momentum conservation 
reads 

k = k 2 + k 3 , (2.55) 
and, in the partonic CM frame, one has 

k 2 = zk (2.56) 
where z = 1 — £,s/(s — M rec ). A factorized expression holds in this case too 



'\ q > 



47ra s 
y^i k 2 ■ k 3 



P qq (z)B qq , =Q 



1 



167TQq 



(l + z 2 ) B 



n 



(2.57) 



The finite quantity needed in the application of the subtraction method is now £ 2 (1 — y) Tl? qq i, 
that is given by 



= c u '~ , u 1 ; 2 

y=l F ZS 



(1 + z 2 ) B qq , . 



(2.58) 



We note again that, in this collinear limit, the contributions [£ 2 (1 — v)T^? qq f] vanish, because 
of the behaviour of the S functions. 

The contribution !Z qq > is also singular when the outgoing gluon becomes soft, i.e. when 
k^ —7- 0. In both the two phase-space parametrizations (^'ad an< ^ ^rad)> this nm it * s 
approached when £ — > 0. The Born process has more than 3 colored particles, so that, 
in general, one may expect that soft singularities factorize in terms of the color ordered 
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Born amplitudes [24]. However, in this case, the color algebra simplifies, because of the 
exchange of an intermediate colorless particle, and we have complete factorization on the 
Born squared amplitude. The 1Z qq i contribution in the soft limit (eikonal approximation) 
is given by 

i n il'} £^0 = 8TO sC P { k h \!l e , h , + 7T , t\n* , 73 " 0(U M .u,2 ) B 11' ■ ( 2 - 59 ) 



The radiation variable y assumes different meaning in the case of ISR or FSR (see sec. 2.2.2). 
In the ISR case, we have the finite contributions 



16 



= 47ra s C F <. 

«=o \s(l±y) (h ■ k 3 )(k 2 ■ h) (h-k 3 ) 2 



where k\ = k 3 /£ identifies the direction of the soft gluon. In the FSR case we have instead 

s{l-y) 4(s-M 2 ) M 2 (l-y)^ 



e(i-y)^ qq 2 , 



= 4ira s C F { s h 



(2.60) 
nsteac 

S 3 ' 2 B qq/ , 



(2-61) 

with £ 2 = Zk^/i/s, defined in the partonic CM frame. 

The i-channel processes Ttbq and 1Z q b are dealt in an analogous way, either for the 
collinear and the soft limits. All the other processes have only ISR collinear singularities: 
the corresponding limits can be obtained from eq. (2.53), substituting the appropriate AP 
splitting kernel and the Born term. 



3. The POWHEG implementation 
3.1 Generation of the Born variables 

In the POWHEG method, we first generate the Born kinematics according to the B function, 
which is the integral of the full NLO cross section at a given value of the underlying Born 
kinematics. It is defined as follows: 



B = £( s ) + -B(t) , 



where 



with 



B 



(s) 



(3.1) 
(3.2) 



qq' 



B qq , (# 2 ) = B qq i (* 2 ) + V qq , (* 2 ) + J d^i Off (* 2> 



+ 



[ d*£3 E (** *55) + R 9 q,(s) (*2, *S3) + R q9 ,(s) (*2, *£5) 

J L © 

^ (* 2 ,e) + Gf ($ 2 , 



1 ^£ 

'x e z 



'Gf (*2,e) + G| 9 (*2,e) 



(3.3) 
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and where 



with 



B (t) = Yl [ B i b + Bb i\ ' 



(3.4) 



B qb (* 2 ) = B qb (* 2 ) + y gb (* 2 ) + | d<&™ R i% 2 (* 2 , 



+ / d*S5 



E^f (* 2 ,<0+^W) (*2,*E3) 



@ 



+ 



»6 

1 r 



G* (* 2>e ) + G£ (# 



2,63, 



+ 



1 cfe 



Gf> (* 2>e ) + Gf (# 2>e )l ■ (3.5) 



The f?^ contribution can be obtained from eq. (3.5) by simply exchanging all flavour 
indexes and substituting ® -H- e. 

According to the POWHEG notation, in eqs. (3.3) and (3.5) we have traded the B, V, 1Z 
and Q quantities with the corresponding capital letters, obtained by multiplying them with 
the appropriate luminosity £, defined in terms of the parton distribution functions (PDF) 

C ff {x e ,x e ) = ff(x 9 ,/4) ff,(x e ,nl) . (3.6) 

All the integrals appearing in the above equations are now finite. In fact, following the 
FKS subtraction scheme, the hatted functions 



and 



•tS FSR 



i 



l 



(it z/) e 2 ^ 



'J 



[(1-y) e 2 ^ SR ] 



(3.7) 



(3.8) 



have only integrable divergences when integrated over $J. ad and $rad respectively. ^ Some 
care should still be taken when dealing with the plus distributions. For more details we 
refer to refs. [24] and [32]. 

Following ref. [24], we introduce the B function, defined such that its integral over the 
radiation variables, mapped onto a unit cube (^{C,y,4>} — > |^ad' ^rad' ^md}) ' gi yes 



B 



rf 3 -^rad B . 



(3.9) 



tin our case, for both the s- and t-channels, 

-pe _ / n>3,e 6 n>3,e p3,e £2,e R 3,e R \ 

K-ij - \K qq > , K qg,(.s) > K qb > U bg ' H bg ' H bq > K 19,(t)j 
T^FSR f 53,2 A3, 2 A3. 2 "1 

l<. i:j - {K qql ,K qb ,K bq | . 
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The generation of the Born variables <&2 is performed by using the integrator-unweighter 
program MINT [43] that, after a single integration of the function B over the Born and 
radiation variables, can generate random values for the variables {$2,^rad}, distributed 
according to the weight B ($ 2 ,X rac i). We then keep the *2 generated values only, and 
neglect all the others, which corresponds to integrate over them. At this stage, we also 
need to choose a Born flavour structure (/& in the language of ref. [24]) with a probability 
proportional to its relative weight in the B function (see eqs. (3.2) and (3.4)). The event is 
then further processed, to generate the radiation variables, as illustrated in the following 
section. 

3.2 Generation of the hardest-radiation variables 

Radiation kinematics is generated using the POWHEG Sudakov form factor. For a given 
underlying Born kinematics ($2) and flavour structure (fb), the Sudakov form factor can 
be expressed as 

Afi>(* 2 ,Pr)= II A£(* 2 ,Pt), (3.10) 
"re{a r |/ 6 } 

where one needs to include in the product all the projected real contributions that have, as 
singular limit, the generated underlying Born. In our case, for the s-channel, we can write 

A™'(* 2 ,p T ) = A«£(* 2>Pr ) A|l' R (* 2 ,p T ) , (3.11) 

where 



A^r(*2,Pt) = exp I - J 



ISR E @ Kf (*s) + Rg* M (*s) + R q9 ,(s) (*s) 



^rad 



x^ t , isr (<I>3)-Pt)> (3.12) 



and 



A|!' r (* 2 ,Pt) = exp|-| d$ F r S J R B 99 q ^ 0(k T , FSR (& 3 ) -p T )| . (3.13) 

For clarity, here we indicate with Rg q '( s ) the real contribution of gq type that corresponds 
to the underlying Born qq' . The functions fc T)ISR ($3) and A; TjFSH ($3) measure the hardness 
of the radiation in the real event. In case of ISR singular processes, we chose as hardness 
variable the exact transverse momentum of the emitted parton with respect to the beam 
axis. In terms of ^ad' * ms ^ s gi ven D Y 

kl^= s l e(i-y 2 ) = ^ ri) e(i-y 2 ). (3.14) 

For the FSR singular processes, instead, we use as hardness variable the exact transverse 
momentum of the FKS parton with respect to the other light outgoing parton, evaluated 
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in the center-of-mass frame. In terms of 3>fad ' ^ ms ^ s gi ven by* 

k 2 T , FS n = ^ 2 (1 " V 2 ) • (3-16) 

The generation of the hardest radiation is performed individually for A 9 g R and Apg R , and 
the highest generated k T is retained. This corresponds to generate according to eq. (3.11), 
as shown in Appendix B of ref. [24]. If k T is below a given cut, p™ m , no radiation is 
generated, and a Born event is returned. 

The upper bounding functions for the application of the veto method have been chosen 
in the following way: 

E @ Kf (*8) + V,W (*») + RwM (*>) j lS n (& ^ < NlS n Mg^ga) . 

J rad (*2, < ±'rad) < ^ an' TTi 77K (A 18 ) 



for ISR, and 



B q< ,(* 2 ) radV z ' rad; " C (1 - 2/ : 



iV(* 2 ) radV z ' rad; " 99 e(i-y 2 ) 

for FSR. 

The same procedures holds also for the t-channel case, with appropriate modifications 
in formulae (3.11)-(3.18). 

The method used to generate radiation events according to these upper bounding 
functions is analogous to the one described in Appendix D of ref. [25], and we do not 
repeat it here. 

As a final remark, we also point out that single-top s- and t-channel Born cross sec- 
tions vanish at some points in the Born phase space, as one can argue by looking at 
eqs. (2.43)-(2.46). For this reason, special care has to be taken during the radiation gen- 
eration procedure. We handled this problem using the same method described in sec. 3.3 
of ref. [29]. We thus refer to that paper for further details. 

3.3 Top- quark decay 

The calculation we have described so far leads to the generation of events with an undecayed 
top quark. We include the decay kinematics effects in an approximate way, by requiring 
that the decay products are distributed with a probability proportional to the tree-level 
cross section for the full production and decay process. This procedure was first suggested 
in ref. [37]. In the following we describe our implementation, focusing upon the decay 
t -»• bW + -»• blv. 

We first generate a Born-like or real-like event according to the POWHEG method. In 
both cases we denote the set of variables that parametrize the undecayed momenta as $ PO w 



* Since for y — > —1 no singularities arise in the FSR case, another possible choice for &t,fsr would be 

fcT,FSR = |C 2 (1"2/), (3-15) 

that has the same behaviour of eq. (3.16) in the collinear limit but has a simpler functional form. We have 
checked that no sizable differences arise if one uses eq. (3.15) instead of eq. (3.16). 
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and the corresponding flavour structure as /. As described at the end of sec. 2.2.1, at this 
stage the top virtuality M 2 is distributed according to a Breit-Wigner function. We write 
the tree-level cross section for production and decay in the following form 



^dec = Ys £ " M dcc( $ POW, <$>t^bt v ) d$dec , (3.19) 



where C is the luminosity factor and M. dec is the squared amplitude corresponding to the 
full decayed process that originates from the undecayed process f} For consistency, the 
squared amplitude M dcc must include only resonant graphs (i.e. graphs where the top 
momentum equals the sum of the b, £ and v momenta). We write the full phase space, 
including the decay, in the factorized form 

d$dec = d$pow d$ t ^ biv , (3.20) 

where $ PO w is the undecayed (POWHEG) phase space and <& t ^biv is defined implicitly by this 
equation. We notice that 

-Kndcc x BR(t ->• blv) = J M f dcc d^ ilv , (3.21) 

where -M undcc is the undecayed squared amplitude, i.e. the Born or real amplitude that we 
used throughout the computation. Thus, the differential probability dP(& t ^ b i„\& POW ) f° r 
the generation of $ t ^biv from a given undecayed kinematics $ PO w is 

dP(* t ^|*Pow) = 77 1 .,- , M ^f™'*»f d * M . (3.22) 
BR(t -> Mv) M{ Qdcc (<S> POW ) 

To generate efficiently <& t -*bi v distributed according to (3.22) we use the hit-and-miss tech- 
nique and so we need to find an upper bounding function for dP. This bound can be 
guessed from the structure of the top decay. In our case, we use as upper bound for the 
ratio A^dc C ($Pow,^t^b&)/-^u n dcc(^Pow), the expression 

Udec\M , = -/Vdec 7775 2\9 , 2772 7772 2 \9 , 2 tW ' (3.23) 

(M 2 - mj) 2 + m\ Tf {Mf v - m 2 w ) 2 + m 2 w Y 2 W 

where M 2 v = (kg + k u ) 2 and Mt^bW and M. w ^i u are the decay squared amplitudes 
corresponding to the subprocesses in their subscripts. In the previous formula, as well 
as in M dec , finite-width effects have been fully taken into account. One can predict the 
appropriate value for the normalization factor iVdec as explained in ref. [37] or compute it 
by sampling the decay phase space <& t ^btv and comparing Ud cc with the exact expression, 
in such a way that the inequality 

■ML^pow, St-*,,) < -^Lcc(^pow) U dcc (M 2 , ^ blv ) (3.24) 
holds. The veto algorithm is then applied: 



§ The full tree-level squared amplitudes M dec have been obtained using MadGraph. 
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1. First one generates a point in the phase space <& t ^btw 

2. Then a random number r in the range [0, Ud ec (M 2 , Q t ^uv)\ is generated. 

f f 

3. If r < M dcc (® POW , ®t^bi v ) /M J nndec (<S>Pow), keep the decay kinematics and generate 
the event. Otherwise go back to step 1. 



4. Results 



In this section we present our results and comparisons with the fixed order (next-to-leading) 
calculation and with the MCSNLO 3.3 and PYTHIA 6.4.21 Shower Monte Carlo (SMC) pro- 
grams.^ We have used the CTEQ6M [44] set for the parton distribution functions and the 
associated value of = 0.226 GeV. Furthermore, as discussed in refs. [24, 25], we use 
a rescaled value A MC = 1.569 A^ in the expression for q s appearing in the Sudakov form 
factors, in order to achieve next-to-leading logarithmic accuracy. 

Although the matrix-element calculation has been performed in the massless-quark 
limit (except, of course, for the top quark), the lower cutoff in the generation of the 
radiation has been fixed according to the mass of the emitting quark. The lower bound 
on the transverse momentum for the emission off a massless emitter (u, d, s) has been 
set to the value p™ m = y/5 A MC . We instead choose p™ m equal to m c or m& when the 
gluon is emitted by a charm or a bottom quark, respectively. We set m c = 1.55 GeV and 
m b = 4.95 GeV. 

The renormalization and factorization scales have been taken equal to the radiated 
transverse momentum during the generation of radiation (see eqs. (3.14) and (3.16)), as 
the POWHEG method requires. We have also taken into account properly the heavy-flavour 
thresholds in the running of a s and in the PDF's, by changing the number of active 
flavours when the renormalization or factorization scales cross a mass threshold. In the 
B calculation, instead, /j,r and \xp have been chosen equal to the top-quark mass, whose 
value has been fixed to nit = 175 GeV. In all the comparisons, we have kept the top-quark 
virtuality M 2 fixed to mf, so that matrix elements have been evaluated assuming F t = 0. 
We have also set Ty/ = in all the propagators. The other relevant parameters are 

M w = 80.4 GeV , sin 2 9§ = 0.23113 , a^(mt) = 127.011989 . (4.1) 

From the above values, the weak coupling has been computed as g = yJAita cra / sin 0$. In 
addition, for sake of comparison, we fixed the CKM matrix elements equal to 



u 

^CKM = C 
t 



d s b 
I 0.9740 0.2225 0.0000 ' 
0.2225 0.9740 0.0000 
\ 0.0000 0.0000 1.0000 



(4.2) 



^This newest update of PYTHIA yields more consistent results when multiple interactions are turned on 
in user-initiated processes (see the release notes in http://projects.hepforge.org/pythia6/). 
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In order to minimize effects due to differences in the shower and hadronization aigorithms, 
we have interfaced POWHEG with the HERWIG angular-ordered shower when comparing with 
MCONLO and with the p T -ordered PYTHIA shower when comparisons with PYTHIA have been 
carried out. 

All the following results have been obtained assuming that the top decays semilepton- 
ically (t — > b£v), as explained in sec. 3.3, but removing the branching ratio, so that plots 
are normalized to the total cross section. 

We present a few distributions, done mainly for comparison with MCONLO and with the 
NLO calculation. Some of them are "unphysical" , i.e., for example, when talking of the top- 
quark momentum p t , we refer to the exact p t taken directly from the MC shower history, 
right before the top decay. For sake of simplicity, we also force the lightest 6-flavoured 
hadrons to be stable after the hadronization stage of SMC programs. 

Jets have been defined according to the k T algorithm [45], as implemented in the 
FAST JET package [46], setting R = 1 and imposing a lower 10 GeV cut on jet transverse 
momenta. We call "top jet" the jet that contains the hardest 6-flavoured hadronjl which 
will, most of the time, come from the top-quark decay. The other reconstructed jets will 
come from the shower of massless partons, and we call them "light jets".** In this way, 
the momentum p t of the top quark and the momentum of the top jet are different, since 
the last may or may not include all the particles from the top decay and shower. 

4.1 Tevatron results 

We start comparing various kinematical variables for single-top s-channel production at 
the Tevatron pp collider. In fig. 2 we have collected the following distributions: 

• In panels (a) and (b) we show the transverse momentum p\ and the pseudorapidity 
77* of the top quark and in panel (c) we show the hardest jet transverse momentum 
p^ 1 . The agreement with the fixed-order calculation and with the MC@NL0 results is 
very good. Only the top transverse-momentum distribution shows a tiny mismatch, 
our result being slightly softer than the NLO and the MC@NL0 ones. When interfacing 
POWHEG with PYTHIA, we instead find full overlapping with the NLO result. It is thus 
likely that this small feature may be attributed to shower effects. 

• In panel (d), we plot p™ 1 ^ 1 , the relative transverse momentum of all the particles 
clustered inside the hardest jet. This is defined as follows: 

— We perform a longitudinal boost to a frame where the hardest-jet rapidity is 
zero. 

— In this frame, we compute the quantity 



"Here we mean precisely b-flavoured, i.e. not fe-flavoured, that arises in the production process. 
**In the fixed-order calculation, instead, the top quark is not decayed, and the top jet corresponds to the 
jet that contains the top quark. 




(4.3) 
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where k l, s are the momenta of the particles that belong to the hardest jet that, 
in this frame, has momentum p ' n . 

This quantity is thus the sum of the absolute values of the transverse momenta, 
taken with respect to the jet axis, of the particles inside the hardest jet, in the 
frame specified above. The plot shows a marked disagreement between fixed order 
calculation and showered results. This disagreement is well understood, since the 
observable we are considering is a measure of the spreading of the hardest jet. Thus, 
its shape is strongly affected by the Sudakov form factor and it is well described by 
SMC programs. The NLO calculation cannot give, instead, a reliable estimate, since 
when Pt 1 ' -7 ' 1 — > the differential cross section diverges. 

• In plots (e) and (/), the next-to-hardest jet transverse momentum p^ 2 , and the 
transverse momentum of the system made by the top quark and the hardest jet, 
Pt^\ are shown. We see a remarkable good agreement between our program and 
MCONLO, while sensible differences with respect to the NLO results are present. At 
the NLO parton level, pi 2 and balance against each other, so that the two 
distributions coincide down to the minimum p T cut present in the first plot. 

In plot (e), we see an enhancement of the showered results at intermediate values of 
p T , while in plot (/) we see a low-p T suppression and an enhancement at interme- 
diate and high p T . The low-p T suppression is clearly a Sudakov effect. The high-p T 
enhancement comes instead from events in which the hardest parton is well balanced 
against the top quark, but where many hadrons, coming from the hardest parton, end 
up in the top jet, and are thus removed, or they end up out of the jet cluster. This 
creates an artificial imbalance, and thus an effective p T for the (tji) system. These 
effects are so pronounced because the cross section for a balanced top-quark-hardest- 
jet system is much higher, since it does not require the production of an additional 
hard parton. We have verified this hypothesis by analyzing POWHEG outputs before 
the showering stage, either clustering or not the b quark coming from the top decay. 
In the case where the b quark is included in the analysis (and the jet containing the b 
is removed from the jet sample), we see a marked rise of the p T tail. A further rise is 
observed when the shower is turned on, and may be attributed to energy lost out of 
the hardest jet cluster due to showering. We see no such effect for the next-to-hardest 
jet spectrum in plot (e). There, the raise at medium p T may be attributed to the 
shower p T smearing. 

• Finally, in plots (g) and (h), the pseudorapidity of the top-quark-hardest-jet 
system and the azimuthal difference A(f) t .j 1 = \(f> t — are shown. The pseudorapid- 
ity of the {tji) system shows an expected discrepancy between the showered results 
and the fixed order one: radiation near the beam axis is suppressed by the Sudakov 
form factor but not in the NLO result, giving rise to the higher tails at large |7/* J1 )|. 
In plot (h), MCSNLO and POWHEG differ instead from the fixed order result for a kine- 
matical reason: at the parton level, having at most three particles, there is no phase 
space for the next-to-hardest jet to recoil against the (tji) system when Acpt-j-^ < n/2. 
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Figure 3: Comparisons between POWHEG, MCONLO and NLO results for t-channel top production at 
the Tevatron pp collider. 
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Figure 4: Comparisons between POWHEG and PYTHIA results for t-channel top production at the 
Tevatron pp collider. 
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A similar set of comparisons is presented in fig. 3 for the i-channel production mech- 
anism, always at the Tevatron. The agreement between POWHEG and MC@NL0 is as good as 
before for inclusive quantities, or even better. In particular, the slight mismatch in the top 
transverse-momentum distribution completely disappears, as one can see in plot (a). For 
all the other plots, considerations similar to the s-channel case remain valid. 

In fig. 4 the same set of plots are shown, comparing POWHEG and PYTHIA. We have good 
agreement for most distributions, after applying an appropriate K factor to the PYTHIA 
results. Only minor differences are present in the high-p T tail of distributions in panels (e) 
and (/). 

As a final comparison, in the left panel of fig. 5, we show , the transverse-momentum 
spectrum of the hardest 5-flavoured hadron, after imposing the rapidity cut |yg| < 3. In 
the i-channel, this hadron will come most probably from an initial-state gluon undergoing 
a bb splitting. The b quark is then turned into a t while the b quark is showered and 
hadronized. We see that, while POWHEG and MCONLO are in a fair agreement in the medium- 
and high-p T range, sizable differences are present at low p T . These discrepancies are most 
probably due to the disagreement that one can notice in the yg distribution (right panel 
of fig. 5), and to a smaller extent to a different implementation of the inclusion of 6-mass 
effects by both programs (just before the showering stage). 




Pt [GeV] y- B 



Figure 5: Comparisons between POWHEG and MCONLO results for the hardest ^-flavoured hadron 
transverse momentum (left) and rapidity (right), for i-channel top production at the Tevatron pp 
collider. Rapidity cuts are highlighted. 

We also plot in fig. 6 the same quantities comparing POWHEG interfaced to PYTHIA with 
respect to PYTHIA alone. A large mismatch in the high-p^ spectrum is clearly visible in 
the left panel. This observable is particularly sensitive to real matrix-element effects, not 
present in PYTHIA. Concerning the low-p^? behaviour, we see that here the difference is 
much less pronounced than in fig. 5. Furthermore, the aforementioned mismatch in the yg 
distribution is no longer present, as one can see in the right panel. 

By comparing figs. 5 and 6, one immediately notices the different behaviours of the 
two Monte Carlo programs that we are interfacing to. We observe that the HERWIG shower 
and hadronization create an enhancement at large values of \yg\, which is not present in 
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Figure 6: Comparisons between POWHEG and PYTHIA results for the hardest 6-navoured hadron 
transverse momentum (left) and rapidity (right), for i-channel top production at the Tevatron pp 
collider. Rapidity cuts are highlighted. 

PYTHIA. This feature is known to the HERWIG authors, ^ and is traced back to a mismatch of 
the scale at which backward evolution is switched off, with the scale at which the 6-quark 
density is turned on in the pdf's. The effect is more pronounced in MCONLO, probably due 
to the fact that POWHEG does not rely on HERWIG for the generation of the hardest splitting. 

4.2 LHC results 

In figs. 7 and 8 similar results are reported for the LHC pp collider. Only plots for the 
i-channel production are shown, the s-channel process having a negligible impact at the 
LHC. 

Figure 7 contains comparisons between POWHEG, MCSNLO and NLO results. No signifi- 
cant differences with respect to what we observed at the Tevatron arise in any plot, so that 
we refer to the previous section for comments. 

In the PYTHIA and POWHEG comparisons shown in fig. 8, we immediately notice that 
the POWHEG enhancement of high-p T tails in panels (e) and (/) is here more marked, even if 
still small. This may again be related to the lack of matrix-element corrections in PYTHIA, 
resulting in larger discrepancies at the LHC with respect to the Tevatron case. 

In panels (c) and (e), one can also notice different low-p T shapes with respect to 
the same plots showing the POWHEG+HERWIG results of fig. 7. We have verified that these 
differences are due to the inclusion of multiple interactions (MI) in the default PYTHIA. $ 
If we limit ourselves to the results without MI (i.e. setting MSTP(81)=0 in PYTHIA), the 
agreement is much better. 



"See M. Seymour's talk in http://bwhcphysics.lbl.gov/vplusjets.html. 

H These account for events where more than one parton pair in the same incoming hadrons give rise to 
hard interactions. 
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Figure 8: Comparisons between POWHEG and PYTHIA results for t-channel top production at the 
LHC pp collider. 
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4.3 Top-quark decay 

As explained in sec. 3.3, in our calculation we have implemented spin correlations in top 
decay. Sizable effects are thus visible when comparing our results with SMC programs that 
do not implement them. MC@NL0 accounts for these effects with approximately the same 
method that we use. Hence, we expect to have good agreement with MCONLO and visible 
discrepancies when comparing with PYTHIA. 

Due to the V-A structure of the weak current, the best observables to highlight eventual 
discrepancies are those involving the angle between the charged lepton t coming from top 
decay and the direction of the down-type quark entering the W vertex involved in top 
production, as shown in fig. 9. 




(a) s-channel (b) t-channel 

Figure 9: Lepton (1) and down-type quark (2) used to study spin correlations in top decay. 

At the Born level, the down-type quark direction is possibly identified with the beam 
axis for s-channel production, while, for t-channel production, it often corresponds to the 
hardest jet axis (see ref. [47] for further details). 

For sake of comparison, we have set the top virtuality M 2 = m 2 and we have taken 
the values Tt = 1.7 GeV and Tyy = 2.141 GeV in the evaluation of upper bounds of the 
decay amplitudes in eq. (3.23) and in the decayed matrix element . Furthermore, we 
have applied cuts similar to those used in ref. [37], both for the Tevatron and for the LHC, 
namely 

p^>20GeV, |r/ B |<2, (4.4) 

pi > 10 GeV, |t/|<2.5, (4.5) 

p v T > 20 GeV . (4.6) 

We denote with the superscript B the top jet, i.e. the jet that contains the hardest 
6-flavoured hadron (not the b). In single-top processes, this comes almost exclusively from 
the bottom quark emerging from top decay. In i-channel production, in order to isolate a 
central hardest light jet, we apply the further cuts 

p^>20GeV, \r] jl \<2.5. (4.7) 

In fig. 10 we show comparisons for the Tevatron pp collider. On the left panel, we plot 
the s-channel differential cross section as a function of cosx, where \ ls the angle between 
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Figure 10: Comparisons between POWHEG, MC@NL0 and PYTHIA angular correlations for s- (left) 
and i-channel (right) top production at the Tevatron pp collider. 



the hardest charged lepton £, which we assume coming from top decay, and the direction 
of the incoming parton with negative rapidity (the e direction of the z axis) , as seen in the 
top rest frame. Such angle is sensitive to the spin correlation between £ and the incoming 
d quark, which, at the Tevatron, is pulled out mostly from the antiproton traveling in the 
negative direction. On the right panel, we plot the i-channel differential cross section as a 
function of cos 9, where 8 is the angle between £ and the hardest jet, always evaluated in 
the top rest frame. In both plots, we observe a remarkable good agreement with MC@NL0 
and the expected discrepancy with PYTHIA, that only performs a spin-averaged top decay. 
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Figure 11: Comparisons between POWHEG, MC0NL0 and PYTHIA transverse momentum and pseudo- 
rapidity of the lepton coming from the top decay, for s-channel top production at the Tevatron pp 
collider. 



In fig. 11 we plot the transverse momentum and pseudorapidity of the hardest charged 
lepton, for s-channel production at Tevatron. The difference between PYTHIA and POWHEG 
(or MCSNLO) can be shown to arise because of spin-correlation effects. To test this, we run 
POWHEG with an undecayed top in the final state, leaving PYTHIA to perform the decay: 
after rescaling the plots with the appropriate K factor, we obtain the same behaviour as 
PYTHIA standalone. 

In fig. 12, the same distributions of fig. 10 are shown for the LHC collider. The same 
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considerations done for the Tevatron apply for the LHC results. 
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Figure 12: Comparisons between P0WHEG, MC0NL0 and PYTHIA angular correlations for s- (left) 
and i-channel (right) top production at the LHC pp collider. 



4.4 Dips in the rapidity distributions 

In previous works [25, 26, 29, 32], we have extensively discussed the presence of sizable 
mismatches between POWHEG and MC@NL0 results in the rapidity difference between the 
hardest jet and the heavy system recoiling against it. More specifically, the MC@NL0 results 
exhibit, for this quantity, a dip at zero rapidity, not visible in POWHEG. This problem was 
originally pointed out in ref. [48] in the framework of tt production, and its origin was 
traced back to HERWIG, that shows an even deeper dip for the same quantity. In ref. [33], 
in the framework of Higgs production, this problem and its Shower Monte Carlo origin was 
accurately studied. 

In single-top production, the suitable quantity where to observe this mismatch is the 
rapidity difference between the top-quark-hardest-jet system and the next-to-hardest jet. 
As one can see in fig. 13, in this case a dip in the central rapidity region is already present 
at the next-to-leading order. This feature may mask an eventual dip in MC@NL0. In fact, 
the two showered results are fairly similar, with the dip being slightly more pronounced in 
MC@NL0. 

In recent talks [49, 50, 51], one of us proposed a possible explanation of the presence of 
these dips in the MCONLO results. In the following we illustrate this explanation and show 
that it is also compatible with the case at hand. 

We can schematically represent the MCONLO cross section for the hardest emission with 
the following formula 



da = B MC ($ n )d$ n 



+ 



pMC/* \ 



A MC (* n ,t ) + A MC ($ n ,t) 



B(*n) 



event 



MC shower 



ii(* n+1 )- J R MC (* n+1 ) 



(4- 



event 
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Figure 13: Comparison between POWHEG, MCONLO and NLO results for the rapidity difference 
between the rapidity of the top-quark-hardest-jet system and the rapidity of the next-to-hardest 
jet, for t-channel top production at the LHC pp collider. Plots are normalized to the total cross 
section. 



The terminology "§" and "H events" is denned in the original MCONLO papers [19, 20]. We 
have 



(4.9) 
(4.10) 
(4.11) 



B MC (3> n ) = S(* n ) + K(#„) + J R MC (& n+ i) , 

i? MC (* n+1 ) = J B(# n )^ip(z), 



A MC (* n ,t) = exp 



2vr t 



where P{z) are the Altarelli-Parisi splitting kernels and (i$ r ad = ^Jid = dz dt deft / (2it) . 
Notice that, on the right hand side of eq. (4.9), divergent quantities appear, and only their 
sum is finite. In the MC@NL0 framework, they are dealt with the subtraction method. 

The "MC shower" factor in eq. (4.8) shows that the hardest emission is produced 
by running the HERWIG shower Monte Carlo, starting with the event kinematics In 
fact, the Monte Carlo may not generate the hardest radiation as its first emission. It was 
shown in ref. [23], however, that formula (4.8) does correctly represent the hardest emission 
probability up to subleading effects, that we here assume to be irrelevant for our argument. 

In the production of a high-p T parton, formula (4.8) yields 



da « 5 MC (* n ) 



R MC (*n 



+1/ 



d* n d*^+ i?(*„+i)- R MC (& 



n+l, 



m , ) ^;:: d + ( - 1 ) i? Mc (* n+1 ) d* n d^ d , 



(4.12) 



0(a s ) 
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where we have used the fact that A MC ($ n ,i) ~ 1 in this limit. The first term correctly 
describes the hard radiation in the whole phase space. The second term, while formally 
subleading in a s , is responsible for the dip. In fact, the dip present in HERWIG propagates 
here with a weight proportional to (B MC /B — 1). Although subleading, this term can be 
significant for processes with large K factors. 

In the processes studied so far, this ratio was significantly higher than 1 (see, for 
example, gg — > H), so that the effect was particularly visible. In single-top production, 
instead, due to the small NLO K factor, one has B MC / B « 1. This, together with the fact 
that the fixed NLO result already presents a central dip for — y-? 2 , results in small 

discrepancies between MCONLO and POWHEG (see fig. 13). 

We notice that a similar mechanism (i.e. via a large B/B factor) for generating large 
NNLO terms operates also in POWHEG, and has been discussed in ref. [32] in the framework 
of Higgs production, as being responsible for a hard Higgs boson p T spectrum. In POWHEG, 
however, this mechanism cannot generate any dip, since here HERWIG has no role in the 
generation of the hardest radiation. 



5. Conclusions 

In this paper we have described a complete implementation of s- and ^-channel single-top 
production at next-to-leading order in QCD, in the POWHEG framework. This is the first 
POWHEG implementation of a process where both initial- and final-state radiation is present. 
The calculation for top production has been performed within the Frixione-Kunszt-Signer 
subtraction approach [38, 39], modified according to ref. [24]. We accounted for spin- 
correlation effects in top-quark decay with a method analogous to the one proposed in 
ref. [37]. The results of our work have been extensively compared with the MC@NL0 and 
PYTHIA Shower Monte Carlo programs, together with the fixed next-to-leading order cal- 
culation, both for the Tevatron and for the LHC. 

The MC@NL0 results are in good agreement with POWHEG, also for quantities sensitive to 
angular correlations in top decay. 

The PYTHIA results, normalized to the total NLO cross section, show fair agreement 
with ours for inclusive quantities that do not involve the top-decay products. As expected, 
we have found sizable mismatches with PYTHIA when considering distributions involving 
top-decay products, such as angular-correlation measurements and charged- lepton trans- 
verse momentum and pseudorapidity. We have also found differences between our results 
and the MC@NL0 and PYTHIA ones in the hardest 6-flavoured hadron transverse momentum 
and rapidity. The high-p T mismatch with PYTHIA may be a consequence of the lack of 
matrix-element corrections in the latter, while we attribute the low-p T disagreement with 
MC@NL0 to the sizable difference that we observe in the rapidity distribution. 

The computer code for the POWHEG implementation presented in this paper is available, 
together with the manual, at the site http://moby.mib.infn.it/~nason/POWHEG. 
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